Spin systems dynamics and faults detection in threshold networks 
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We consider an agent on a fixed but arbitrary node of a known threshold network, with the 
task of detecting an unknown missing link/node. We obtain analytic formulas for the probability of 
success, when the agent's tool is the free evolution of a single excitation on an XX spin system paired 
with the network. We completely characterize the parameters allowing for an advantageous solution. 
From the results emerges an optimal (deterministic) algorithm for quantum search, therefore gaining 
a quadratic speed-up with respect to the optimal classical analogue, and in line with well-known 
results in quantum computation. When attempting to detect a faulty node, the chosen setting 
appears to be very fragile and the probability of success too small to be of any direct use. 



I. INTRODUCTION 

Searching and traversing graphs with the use of a quan- 
tum dynamics (discrete or continuous) is a topic of wide 
study. The main findings on the algorithmic side exhibit 
a quadratic gain (producing then a class of Grover-like 
results), when used to search marked vertices in hyper- 
cubes, lattices, and more general objects [HI) and an 
exponential one when transferring information to a dis- 
tant site by free evolution Q . With essentially the same 
formalism, the dynamics of spin systems forms the basis 
of concrete proposals for implementing communication 
buses in a variety of nanotechnology devices 0. The 
central point is always a unitary operator reflecting the 
topology of a graph whose vertices are encoded in pure 
states or a (time-independent) Hamiltonian describing 
the interactions of particles in the actual physical net- 
work, and thus representing a noiseless quantum channel. 

Considering these two contexts together, we take a par- 
simonious agent located on a specific node of a known 
threshold network, attempting to detect an unknown 
missing link/node. We obtain analytic formulas for the 
probability of success, when the missing link/node is 
searched by letting evolve a single excitation from the 
given vertex. For specific parameters, we have an op- 
timal algorithm for quantum search in a deterministic 
fashion. Our working ground is an XX (or, equivalently, 
isotropic XY) system with homogeneous couplings and 
site dependent magnetic fields [3[ (a proposed setting for 
distributed implementation consists of cavities connected 
by optical fibers Q). 

Threshold graphs are applied to the synchronization of 
parallel processes, set packing, scheduling, and to define 
the Guttman scale in the area of statistical surveys [Io| . 
Formally, a threshold graph can be constructed from the 
one-vertex graph by repeatedly adding a single vertex 
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of two possible types: an isolated vertex, i.e. a vertex 
without edges; a dominating vertex, i.e. a vertex con- 
nected to all other vertices. The creation sequence of a 
threshold graph G on n vertices is a word x (G) G {0, 1}™ 
recursively defined as follows: x (G) — x(G — i)y, where 
y = if i is isolated and y = 1 if i is dominating. A 
threshold graph is then characterized by its creation se- 
quence, a property with important consequences. In fact, 
the information required to store a threshold graph on 
n vertices is at most n bits, two threshold graphs are 
isomorphic if and only if they have the same degree se- 
quence, and the recognition can be done in linear time. 
Even if the above definition is somehow restrictive, it is 
valuable that one can construct arbitrarily large thresh- 
old networks with approximately any prescribed degree 
distribution, including the scale- free one @. 

The dynamics of the first excitation sector in the XX 
model is governed by the (combinatorial) Laplacian Q. 
The Laplacian of a graph G = (V, E) is the matrix L with 
entries [L]i j = d (i) if i = j, [L]i j = — 1 if {i, j} G E and 
[L]i,j = if {i,j} i E; here d(i) = \{j : {i,j} E E}\ 
is the degree of i G V. Let Ai > A2 > • • ■ > A n = 
be the eigenvalues of L arranged in the nonincreasing 
order. A graph is a threshold graph if and only if 
Xj = \{i : d (i) > j}\, with j = l,...,n [5j. In other 
words, the degree sequence and the nonzero eigenvalues 
of a threshold graph G are (Ferrer's) conjugate partitions 
of 2 \E\. For example, let us consider the word 001 101 1 . 
The graph G has degree sequence {2,4,4,5,5,6,6} and 
eigenvalues {7, 7, 6, 6, 4, 2, 0}. 

Let \tp t ) = U t \ipo}, where U t = exp(-iit), for tel+, 
and |-0o) is an element of an n dimensional Hilbert space, 
with standard basis {\i) : i € V}. Let pa(i,j,t) — 
\{j I Ut \i) I be the probability that a single excitation trav- 
els from node i to node j after a free evolution of duration 
t (this is also called fidelity); we say that there is perfect 
state transfer when pq (i,j,t) = 1, for i =/= j. Since the 
Laplacian spectrum of a threshold graph is integral, it 
follows that the dynamics governed by L is always peri- 
odic [2, 4]. i.e., there is t such that po (t, i, t) = 1. for 
every i. 
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Firstly, we propose in Theorem [2] a complete charac- 
terization of the class of threshold graphs allowing for 
perfect state transfer. Secondly, with a direct applica- 
tion of this result, we describe how the evolution can 
be used to detect a missing link. The procedure gives a 
straightforward algorithm for (optimal) quantum search 
in certain threshold networks. Indeed, as a special case, 
we have a deterministic version of Grover's search 
Theorem [6] shows that the result can not be extended 
to a large subfamily of threshold graphs. As expected, 
the graphs for which the result holds can be associated to 
modulated chains of length 3, where the couplings are de- 
termined by the number of edges connecting vertices from 
certain subsets (namely, an equitable partition). This is 
a common situation when studying perfect state trans- 
fer 0. Interference effects drive the evolution in a way 
that does not seem to be exploitable for searching (by 
using out technique), unless the missing links belong to 
the subsets. 



Theorem 2 Let G be a threshold graph. When 
G = T(m 1 ,...,m 2 k) (resp. G = T(n%i, . . . , m 2 k+i)), 
PG {i,j,t) = 1 if and only if '{i, j) = (1,2) and in addition: 
t = 7p mi = 2; ni2 = 2 mod 4; and mj = 0mod4, j — 
3, ...,2k (resp. j = 3, . . . , 2k + I). 



Proof. Suppose that we have mi, . . . ,m 2 k £ N with 
mi > 2, and consider the Laplacian matrix L of 
r(mi, . . . , m-ik)- For each 1 < / < 2k, let 07 = Yl P =i m v 
Fix an odd index j between 1 and 2k, and note that if 
rrij > 2, then for any vector u € W mi such that u _L l mj , 

the vector [Oo-^ |u|0 (T2fc _ (Tj ] T is an eigenvector for L cor- 
responding to the eigenvalue Ao(j) = mj + i+m,j + 3 + . . .+ 
m,2k- Letting u\, . . . ,u mj -\ be an orthonormal basis of 

(lmj)" 1 ) we find that Yd=!^ u l u T = I-^-Jm^mj- Note 
also that if j is odd and 2 < j < 2k, then the vector 



II. FAULTY LINKS 

For any p6N, let K p denote the complete graph on p 
vertices, and let O p denote the empty graph on p vertices. 
Let p and l p denote the zero vector and all ones vector, 
respectively, both of order p. For p, q £ N, we use J pyq to 
denote the p x q all ones matrix; often the subscript will 
be suppressed when the order is evident from the context. 
We need two operations: the union and the joint. The 
union, denoted by U, consists of taking two graphs and 
looking at them as a single one whose connected compo- 
nents are exactly the graphs. The join, denoted by V, is 
the graph obtained by taking the union of two graphs, 
plus edges connecting all their respective vertices. The 
following lemma can be deduced from the basic proper- 
ties of threshold graphs and it does not require a proof 
(see, e.g., Q)- It is our main technical tool. 

Lemma 1 Let G be a connected graph on at least two 
vertices. Then G is a threshold graph if and only if one 
of the following two conditions is satisfied: 

(i) there are indices mi, . . . ,m 2 k G N with m\ > 2 
such that G=(( ( (O mi V K m2 ) U O m3 ) V K mi ) . . .) V 



Y(m 1 ,...,m 2 k); 



(ii) there are indices mi, . . . , m2k+i € N with m\ > 2 
such that G = {{{{K mi U m2 ) V K m3 ) U O mi ) . . . ) V 
K rn 2k+1 = r(mi, . . . ,m 2 k+l)- 

The vertices 1., .., m\ correspond to the first subset, 
mi +1) fni +fn2 correspond to the second subset, 
etc. 

The next theorem is the central result of this section. 
It gives a complete characterization of threshold graphs 
with perfect state transfer. As it is usually done when 
studying this topic, we write explicitly the eigensystem of 
the Laplacian. The parameterization is the one described 
in Lemma [1] 
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is also an eigenvector for L corresponding to Xo(j) that 
has length one and is orthogonal to the eigenvectors con- 
structed above. Similarly, if j is an even index between 
2 and 2k, and u\, . . . ,u m .-\ is an orthonormal basis 

of (l mj ) x , then the vectors [0 (Tj _ 1 |u;|0 (72fc _ (Tj ] T , with 
/ = 1, ...rrij — 1, are eigenvectors for L corresponding 
to the eigenvalue Ai (j) = aj + m J+ 2 + nij + 4 + . . . + ni2k- 
Note further that the vector in Eq. ((TJ) is an eigenvector 
for L corresponding to Ai(j). 

Finally, we observe that \ ar>k j ^Jo~2k is a null vector for 
L. It now follows that for each index j between 2 and 
2k, the matrix 
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with x = 0, 1, is an orthogonal idempotent {i.e., a projec- 
tor) for L corresponding to the eigenvalue Ao(j') or \i{j) 
according as j is odd or even, respectively. Here the (1, 1) 
block is of order o-j-i, the (2, 2) block is of order nij, and 
denotes a zero block whose order is determined from 
the context. It now follows that for any real t > 0, we 
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can use the orthogonal idempotents to write 



itAo(l) 


I - 




<mi 














(2) 



E 

j>3,j, odd 



e -«*oU)p[Ao(7)] 



j>2,j, even 



(We note that an equivalent expression for f/ t appears in 
Theorem 3.1 of Q.) 

Now we consider an off-diagonal entry z of Uf, for con- 
creteness, we take z to be in the upper triangle of the 
matrix, and we suppose that z is in a column that lies 
in the jo _ th subset of the partitioning that is induced by 
the indices mi, . . . , rri2k- (We observe in passing that all 
such offdiagonal entries are equal.) From Eq. ©, we 
find that 



ax 



Jo 



2 k 

E ■ 

j=jo+l 



-it\ Pj (j) _ 



(Ti-iaj 



1 

0~2k 



where for each index j, pj is or 1 according as j is odd 
or even, respectively. Consequently, 



z < 



2k 



3=30 + 1 



a,_iCT,' 



1 

0~2k 



Note that m 2k /{a 2 k-io- 2k ) + V^fe = l/aik-i- Thus, the 
summation giving \z\ telescopes, so that 



2k 

E 

3=30 + 1 



cr,_iCT,- 



1 

0~2k 



Hence, \z\ < 2/<jj . In particular, since mi > 2, we see 
that |z| < 1 if jo > 2. 

Suppose now that \z\ = 1. Then necessarily j = lj 
and, since 1 < 2/mi, we find that mi = 2. Indeed, 
it must also be the case that z is the (1, 2) entry of Ut- 
From the consideration of the equality case in the triangle 
inequality, we find that \z\ — 1 if and only if, in addition 
to mi — 2 we have: 

(i) e-^oM = -1, while 

(ii) e -' tA «0) = 1, for all odd j > 3 and 

(iii) e - *** 1 ") = 1, for all even j. 

From (i) we have e -2 * 1 = — 1, while from (ii) it follows 
that e _4 *" l2! = 1 for each I = 2, . . . , k. Applying these 
conditions, in conjunction with (Hi), it then follows that 



e -itm 2 _ _ 1; whne e - it m 2! 

Hence there are integers pi , 



= 1, for I = 1, 
■ , P2k such that 



,k - 1. 



7T 



2pi 



1 



2p 2 



1 



2p,- 



for j = 3, . . . , 2fc. We thus find that m 2 = 
(4p 2 + 2) / (2ft + 1) while ra 3 = %/ (2p t + 1), j = 
3, . . . , 2k. We now deduce that m 2 = 2 mod 4, and 
mj = mod 4 for j = 3, . . . , 2fc. 

Next, we consider the values of t e [0, 27r] for which the 
(1,2) entry of Ut has modulus one, and note that since 
L has integer eigenvalues, an offdiagonal entry of Ut has 
modulus one if and only if the corresponding entry in 
exp(— i(t + 2qir)L) has modulus one for any integer q. 
From the foregoing we find that the only possible values 
of t € [0, 2tt] for which the (1, 2) entry of U t has modulus 
one are t = ir/2 and t = 3ir/2. Since m2 = 2mod4 
and mj = mod 4, for j = 3, ... ,2k, we find readily that 
conditions i)— iii) hold when t = ir/2 or t = 3tt/2, so that 
the (1, 2) entries of exp(— inL/2) and exp(— HttL/2) have 
modulus one. The second part of the theorem follows 
along the same lines. ■ 

Corollary 3 Let G be a threshold graph on n = 4m ver- 
tices with mi — 2 and m 2 = Am — 2. The pair {i, j}, with 
i and j in the set of size mi, can be found with certainty 
by the use ofO(n— 1) evolutions induced by L, each one 
of time 7r/2. 

Proof. Consider K n , the complete graph on n vertices. 
A special case of Theorem [2] is K~ := K n — {i,j}, for 
two distinct arbitrary vertices i and j. When n = Am 
(m > 1), p K - (i,j, tt/2) = 1 and p K - (k,l,n/2) = 0, 
for every pair {k, I}, with k,l ^ Since |£^(i^ n )| ^ 
n 2 for n —> oo, n — 1 applications of U^/2 implement a 
deterministic and optimal version of Grover's search [l[. 



In practice, we are given a network modeled by a com- 
plete K n that is missing a single unknown link {i, j}. By 
starting the dynamics on each possible vertex (or, equiv- 
alcntly, particle) of K~, we can determine the missing 
edge by letting the entire system evolve for a time it/ 2 
and then perform a projective measurement at the same 
vertex. The dynamics governs the behaviour of a walker 
on the network. If the walker has moved, then the new 
position is vertex j and the link {i,j} is missing. 

The same reasoning may be applied to find a miss- 
ing matching. Recall that a matching is a set of vertex- 
disjoint edges. A matching is perfect if it includes all 
vertices. 

Corollary 4 Let K n be a graph on n = Am vertices. A 
deleted matching from K n of size k < 2m can be found 
with certainty by the use of exactly n/2 — 1 evolutions 
induced by L, each one of time tt/2. 

It is just matter of running the dynamics from an ar- 
bitrary vertex i and detecting the missing edge {i, j} in 
the matching. Then, we pass to a vertex k ^ i, j, and so 
on and so forth. A deterministic search requires at most 



m 2 



J2i=3- odd* = n /4 — 1 steps. For threshold graphs, The- 
orem [2] specifies completely the detectable links. From 
the perspective of a direct application to network search 
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based on the described method (free-evolution on a spin 
system), the theorem shows that the complete graph is 
the only threshold graph in which every deleted link can 
be actually found. For all other threshold graphs there 
are some undetectable links. Even if the same proce- 
dure can be certainly generalized to any graph, since 
the search is performed in cliques (i.e., complete sub- 
graphs) , a complete knowledge of cliques is necessary and 
a mechanism to induce evolution only in desired cliques 
is needed. Still, it is useful to remark that implementing 
such a mechanism (which would turn on and off the in- 
teractions between different cliques) permits the transfer 
of an excitation to any desired vertex. 



III. FAULTY NODES 



above, the value cos(7r/a) is attainable for some t. If m 
is odd/even then 

cos(7r/a) = {— cos 2(cn), — cos a(cir), cos(a — 2)(c7r)} , 
{— cos2(&7r), — cosa(&7r), cos(a — 2)(&7r)} 

respectively; in a similar way, if m is even/odd 

cos(7r/a) = {— cos2(c7r), cosa(c7r), — cos(a — 2)(c7r)} , 
{— cos2(&7r), cosa(&7r), — cos(a — 2)(&7r)} . 



Theorem 6 Let mi, . . . , rri2k € N such that mi = 2. 
777 2 = 2 mod 4, and mi = mod 4, for I — 3, ... ,2k. Let 

i + j and g = U + + 



Here we consider the problem of the previous section 
but for vertices. Given a threshold graph G, let G = 
G—j, for some vertex j g V, and let U t = exp(—iL(G)t), 
We ask whether it is plausible that if G has perfect state 
transfer, we are then able to detect the missing vertex. 
The idea is based on taking advantage of the graph struc- 
ture, apart from the edge between the two vertices in the 
set of size m\ . After a technical lemma, we will prove that 
the offdiagonal entries of U t , while bounded away from 1 
in modulus, can have relatively large magnitude, some- 
thing which does not help to give specific information 
about the deleted vertex. As a consequence, we do not 
obtain sufficient information about j. The negative re- 
sult highlights an important role for special symmetries. 
We shall give a proof for G — r(mi, mz, . . . , m2t)- The 
theorem for the case T(m\, 777,2, ■ • ■ , W2fe+i) has a parallel 
statement. 

Lemma 5 Let a £ N be odd. Then 

(i) max{min{— cos 2?;, — cosai, cos(a — 2)f}|t £ [0,27r]} 
= cos(7r/a); 

(ii) max{min{— cos 2?;, cosat, — cos(a — 2)f}|t G [0,27r]} 
= cos(7r/a). 



Proof. Set a = 2m + 1. If 

teS = 



m m + 1 

7T, 7T 



2m + 1 2m + 1 



3m + 1 3m + 2 

7T, 7T 



2m + 1 2m + 1 



then — cos2t < cos(7r/a). On the other hand, if t £ S 
then |cos(a — 2)t\ < cos (n/a). Thus, we find that 

min{— cos 2t, — cos at, cos(a — 2)t}, 

7T 

min{— cos2rj, cosai, — cos(a — 2)t} < cos — , 

a 

for all t € [0, 2?r]. Let b = (m + 1) / (2m + 1) and c = 
mj (2m +1). Next we show that for each of the functions 



(i) If G = T(mi - l,m 2 , . . . ,m 2fc ) then \[U t ]i. \ < 
2/(m 2 + l); 



(ii) If I is even and G — r(mi, 
l,m ;+ i, . . . ,m 2 fe) then 



[U t ]i,j\ < [9 ~ ( 1 ~ cos ; ^ 2k 



,mi-i,mi 



m 2 



0"icr 2 cr 2fc 



1/2 



m 2 



even mi , 
1 



ClCT 2 (T 2 fc 



(3) 



(Hi) If I > 3 is odd and G = T(rn±, . . . ,mi-\,mi 
l,m/ + i, . . . ,m 2fe ) then 



\Pt]i,j\ <[g-\ 1-coB 



m 2 m /+ i 



1/2 







odd m i 


- 1 J 


1 


m 2 


mi+i 


+ 2- 


(Tl(T 2 


o~icri + i 



(4) 



Proof. Let z be an offdiagonal entry of Ut- As in the 
proof of Theorem [21 we find that if z ^ [^t]i,2, [£^2,1 
then \z\ < 2/<i2 < 2/3. Thus, for the remainder of the 
proof, we assume wlog that z is in the (1,2) position. 
Recall that 07 = 2p=i m v 

(i) Note that T(mi — 1, m 2 , . . . , m 2 fe) = T(m 2 + 
1, 7773, . . . , m 2 fe). From the proof of Theorem [2j it follows 
that \z\ < 2/ (m 2 + 1). 

(ii) Suppose that I is even and that G = 
T(mi, . . . ,m ; _i,m ; - l,m ;+ i, . . . ,m 2fc ). Set a = 1 + 
77i 2 + 7774 + . . . + m 2 fe . Suppose that we have posi- 
tive parameters ct/3,7 such that a > (3, 7, and note 



that 



/i 2 



7 



i -ae""'"- 2 ) + Pe~ lta + 7 | = a 
2a7Cos(a — 2)t + 2f3j cos at — 2af3cos2t. From Lemma 
for each t, one of — cos 2t, cos at and — cos(a — 2)t 
is bounded above by cos(7r/a). It now follows that 

|_ ae -«(a-2) + p e -ita + 7 | < ( a 2 + 02. + + 2al + 

2/3 7 cosf +2a/3) 1 /2 = (( a + p + 7 )2 _ (i _ cos ^j) 1 / 2 . 
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Taking a = 1/crj = 1/2, /3 = m 2 /(o-i<7 2 ), and 7 = l/cr 2 fc, 
applying the triangle inequality as in the proof of Theo- 
rem [21 we have 



\z\ < 



t(a— 2) 



< 



i -it(a- 
2 

m 2 



Cicr 2 

™2 _ ita , 



(Tlf7 2 



1 

H 

0"2fe 

1 
C 2 fc 



2k 

y, TTij 
.7=3 J J 



+ 1 



0"lCT 2 



1 

C 2 fe 



which implies Eq. (|3]). 

(m) Suppose that Z > 3 is odd, and that G = 
T(mi,...,mi-i,mi - l,mj+i, . . . ,m 2 fe). Set a = 
?77i + 7713 + ... + m; — 1. As in the proof of 
(ii), we consider positive parameters a > /3, 7, and 
note that |-ae _ltA ° (1) + /3e _ltAl(2) + 7e -' ;tA i( i + 1 ) | = 
|— a + (3e~ 2lt + je~ tat \. As in (ii), we deduce from 
Lemma [5] that |-a + /3e _2lt + 7e _wt | < ((a + /3 + 
7) 2 - (1 - cosf )/?7) 1/2 . Taking a = l/oi = 1/2, 
/3 = 777,2 /(ctict 2 ), and 7 = TO ;+ i/((7icr ;+ i), Eq. (gj fol- 
lows. ■ 

Notice that in the context of Theorems HI the (1,2) 
entry of U t can have large modulus. For example, 
it is straightforward to show that if m,\ = 2, m 2 = 
2 mod 4, and to; = 0mod4, for I — 3, ...,2k and 
G = r(TOi,m 2 ,...,TO2fe_i,TO2fe - 1), then |[J7§]i, 2 | = 

1 /2 

(1-2 ((a 2 k - l)/°2fc)) ■ Similarly, if m 1 = 2, m 2 = 
2 mod 4, and m; = 0mod4, for I = 3,...,2fc + 1 and 
G = r(mi,m 2 ,...,m 2 fe,m 2 fc + i - 1), then |[?7f]i, 2 | = 

{\-2{{a 2k+1 -l)/al k+1 )) 1/2 . 

IV. CONCLUSIONS 

We have shown that the dynamics on a network gov- 
erned by the Laplacian, seen as the Hamiltonian re- 
stricted to the single excitation sector of an XX spin 
system, can detect and find a faulty link in the complete 
graph with a quadratic gain with respect to a determin- 
istic method, and we have extended the observation to 
matchings. The result gives a way to perform optimal 
quantum search (on the complete graph), and a new in- 
sight into the related algorithms. Essentially we have a 
reinterpretation of a fact discovered in Q. Our contri- 
bution is to have put the statement in a more general 



mathematical context, by giving a complete character- 
ized of threshold graphs with perfect state transfer with 
respect to the XX model. Dealing with deleted nodes, 
we have shown that our method does not give any clear 
advantage. The basic idea of the paper is to use a free 
quantum evolution to search a missing item on a net- 
work. The method is different from the ones studies in 
Jll] substantially because the algorithm does not require 
any control. After the set-up of the network, the system 
evolves without intermediate operations. The process is 
distributed because the nodes of the network are identi- 
fied with spin particles. The measurements are local in 
the sense that they are concerned with the single sites, 
independently. That is why we can locate autonomous 
agents on the sites. 

We have considered threshold graphs only, since these 
have integer Laplacian spectrum, and so possess peri- 
odic dynamics, i.e., a necessary condition for perfect 
state transfer (see 0). The vertices involved in the 
phenomenon are special. The symmetry with respect 
to these vertices can be exploited to create a "reference 
point" inside the graph. Looking ahead, one direction 
for further exploration is to determine what information 
about the topology of a spin system (paired with a net- 
work) can be obtained by a free evolution and final local 
measurements. Here, more than designing search algo- 
rithms, it is a matter to determine what kind of graphs 
have some sort of searching capability embedded in their 
structure. The method of the paper can be generalized to 
searching a missing link (or, equivalently, a marked link) 
in any network, with steps involving one clique at a time. 
The method works for vertices if we employ the notion of 
the line graph. However, we have shown that the direct 
detection of a missing vertex is not a natural task for 
the studied dynamics, at least on threshold graphs. It 
is an open question to determine whether the dynamics 
can help to find marked nodes, when taking a different 
initial probability distribution, and if we can obtain the 
quadratic speed-up in this case. 
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